Spectral flow in vortex dynamics of d-wave superconductors 
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We calculate Ohmic and Hall flux flow conductivities in the mixed state of d-wave superconductors 
taking into account the spectral flow along the anomalous branch of the states bound to vortex cores. 
Two parameters governing the strength of the spectral flow at low temperatures are found. They 
are determined by the ratios of the width of the localized levels 1/r and 1) the interlevel spacing 
in the core, and 2) Aq/Ef- These parameters determine three regimes of vortex dynamics. In 
the moderately clean limit the contribution of the spectral flow to the conductivities coincides with 
that for s-superconductors. In the superclean regime in low magnetic flelds both Ohmic and Hall 
conductivities acquire r-independent universal values. The contribution of the spectral flow in the 
superclean regime is suppressed in higher flelds. We also discuss the case of higher temperatures. 



I. INTRODUCTION 



+Ck X (v„ - vl). 



(1) 



X3 

O 
o 



> 

00 

m 
o 
l> 
o 
l> 

ON 



I 

O 

o 



The problem of vortex motion in superconduc- 
tors/superfluids and flux-flow Hall effect attracted atten- 
tion of researchers for years. The renewed interest to this 
problem is due to recent developments in the physics of 
the Hall effect in high-temperature superconductors and 
in the investigations of the superfluid '^He. A microscopic 
theory of mutual friction (friction between the superfluid 
and normal components in presepice of quantized vor- 
tices) was developed long time ago.tJ This theory (see also 
Refs. enables us to calculate the forces on a vortex 
(or, equivalently, Ohmic and Hall conductivities) as func- 
tions of microscopic and external parameters. The tem- 
perature deoendence of the forces on a vortex predicted 
theofeticallya is in good agreement with the experimental 
datal] for ^He-B. The theory also proposes an explana- 
tion for the sign reversal of the Hall angle as an intrinsic 
property of the mctex flow in conventional and high-Tc 
superconductors.aci This theory appHes to the regime of 
free flux flow when pinning is not important. In super- 
conductors this region can be reached experimentally ei- 
ther by preparation of high purity untwinned samples or 
by using high density driving currents. 

However, all the previous calculations were performed 
for the isotropic s-pairing. Since this theory is used for 
analysis of the data for high-Tc materials it is necessary 
to understand whether its arguments are sensitive to the 
symmetry of the pairing state and the presence of gap 
nodes. 

The total non-dissipative force on a vortex can be sub- 
divided into the following three contributions of different 
physical origin (see Ref. H and references therein): 

F = pK X (Vi - Vs) + PnK X (Vs - V„) 



where p, p„ are the densities of the fluid and the normal 
component, respectively, k is the vector of the circula- 
tion of the superfluid velocity around the vortex directed 
along the vortex axis, k = ■nfi/m for a 1-quantum vortex. 
v„ and -Vs are the velocities of the normal and superfluid 
components far from the vortex core and is the veloc- 
ity of the vortex (the vortex lattice). The first term in (|l|) 
is the Magnus force describing the transfer of momentum 
from the superfluid vacuwn to the vortex. The second 
one is the lordanskii forceB due to asymmetric scattering 
of excitations by the vortex core. These first two terms 
could be rewritten as /? x [ps(vs — v/,) -I- p„(v„ — v^)] 
and thought about as contributions of the superfluid and 
normal components. 

In this article we concentrate on the last term in (|^). 
This term describes the transfer of momentum between 
the heat bath and the bouBd. states in the vortex core spe- 
cific for Fermi superfluids|l3 which is due to the spectral 
flow (SF) along the anomalous branch of bound states. El 
For s-wave superconductors the coefficient C in front pf 
this term depends on the ratio of the interlevel distanccEll 
in the core ojq ^ 1^ jEp to the inverse relaxation time 
T for the bound states. When the width of the levels 
l/r is greater than the interlevel spacingjJx^o , the spec- 
trum is continuous and the SF is allowedjj In this limit 
C — Cq, which is the density of the fluid in the nor- 
mal state (Co = mkp/2'!T in 2D, Co = mk^/Sir'^ in 3D). 
The difference p — Cq ^ plA/Ep)"^ is small and de- 
pendence of the total force (|]) is weak. In the opppsite 
limit ujqt 3> 1 the spectrum is discrete and the SFllil is 
suppressed, C — > 0. 

Another explanation of the dependence of the streagth 
of the spectral flow on t was proposed by Stone.Q He 
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showed that quasiclassical bound states inside the vortex 
core rotate with the angular velocity f^nd this rota- 
tion suppresses the spectral flow. However, relaxation 
effectively prevents the rotation if the relaxation time r 
is less than the period of rotation ^-k juj^. Therefore, the 
SF is suppressed in the limit of large loqt and allowed 
in the opposite limit. In the regime of free flux flow the 
only contribution to the dissipative force is also due to 
the bound states. The dissipative force is small in both 
limits and is maximal at wqt ~ 1. 

At finite T ~ Tc a contribution of the continuous spec- 
trum of states above the gap to C is also essential. Since 
the interlevel distance above the gap is given by the cy- 
clotron frequency Wc = |e|^f7^*c ~ bjQ{H / Hc2), this con- 
tribution of the continuouai3 spectrum is suppressedQ in 
the limit ojcT 3> 1. 

In this article we address the question of applicability 
of the theory to the d-wave casejla namely, to (1^2 _y2- 
pairing state, which is believed |-tp represent the order 
parameter in high-Tc materials,ll3 though the analysis 
for an arbitrary pairing state with points of gap nodes 
in 2D (or lines in 3D) would be similar. The structure of 
the states bound to vortex cores in d-superconductors is 
different from that in s-wave materials.LJ In particular, 
low-energy states are extended in directions of gap nodes. 
In addition, at energies below the maximum value of the 
gap Aq bound states coexist with the states in the con- 
tinuous spectrum. So, one can expect new features in the 
contribution of these low-lying levels to vortex dynamics. 

The interlevel spacing in cores of d-wave vorticealj 

is much less than in s-superconductors: lOc sq <ti 
Aq/ Ep if iJ <C Hc2- Comparison of eq to the width 
l/r of the levels suggests the dimensionless parameter 
eQT governing the strength of the SF. 

On the other hand, the angular velocity of rotation of 
Andreev bound states ivoiO) is angle-dependent in d-wave 
superconductors (it follows the angle dependence of the 
gap; see Section 2). Its maximal value w™""^ on the order 
of Aq/Ef (Aq is the maximal value of the gap) is much 
greater than eo- So, this approach proposes a different 
scale for l/r as a boundary between the regimes of ex- 
treme and suppressed SF. In addition, near the gap nodes 
the interlevel distance uJo{9) is verit-small (see Section 2), 
i.e., the density of states is largell3 So the contribution 
of this region needs special attention. 

One of our purposes is to understand, which parameter 
governs the strength of the SF in the d-wave case. In par- 
ticular, this question is important for analysis of exper- 
imental data: one should locate the boundary between 
the regimes of extreme and suppressed SF. In addition, 
it is interesting to estimate the force on a vortex and 
Ohmic and Hall conductivities in various limiting cases 
to compare the results to experiments and previous theo- 
ries. In particular, we look for results different from those 



for s-superconductors. We found new regimes of vortex 
dynamics in the superclean limit in superconductors with 
gap nodes (d-wave), different from those in isotropic (s- 
wave) superconductors. The SF force depends on both 
parameters discussed above, cqt and Aqt/Ef (as well as 
on LOcT at finite T). 

It turns out that it is convenient to explore this prob- 
lem using the kinetic equation approach developed re- 
cently by Stonea for s-wave vortices. We apply it to the 
d-wave case and find the solution of the equation. To use 
the kinetic equation (KE) one should first investigate the 
spectrum of the bound states (Section 2). Then, in Sec- 
tion 3, we solve the KE and discuss its implications on the 
force in various limiting cases in Section 4. We consider 
the 2D case (applicable directly to thin films and pan- 
cake vortices). Generalization to 3D is straightforward 
(see discussion). 

In Section 5 we explore the KE in terms of the density 
matrix in the basis of exact eigenstates of Bogolyubov- 
de Gennes (BdG) hamiltonian, and evaluate parameters 
governing the strength of the SF in a different way. Fi- 
nally, we summarize and discuss our results in Section 
6. 



II. SPECTRUM OF THE BOUND STATES 

In this Section we describe states inside and near vor- 
tices in d-wave superconductors. We consider a ma- 
terial with a circular particle-like Fermi surface (gen- 
eralization to the case of a hole-like Fermi surface or 
coexisting hole- and particle-like parts of a Fermi sur- 
face is straightforward). We assume that the coherence 
length ^0 in the superconducting state is much greater 
than the inverse Fermi momentum kp^ (or, equivalently, 
Aq ^ Ep). Even in high-Tc materials, where Ep /Aq 
is not very large, it is still of the order of 7-10. There- 
fore we can use the quasiclassical approach and solve the 
BdG equations in Andreev approximation. We also as- 
sume that the magnetic field H <C Hc2 so that vortices 
are far apart (and, of course, H > Hd so that the sample 
is in the mixed state). 

In dj.2_y2-wave state the distribution of the order pa- 
rameter in the vicinity of a vortex is given by 

A{r,9)^Aoir)e^^ikl-k!) (3) 

where a and b are the crystal axes, r, ip are polar coor- 
dinates in the position space. We put a = x, b — y. The 
absolute value of the gap Aq (r) equals zero at the center 
of the core and saturates at the bulk value Aq far from the 
core, at distances r ^ ^o- In the bulk A{6) = Aq cos(26'), 
where the angle 9 is defined by k = kp {cos 9; sin 9), and 
there are 4 gap nodes (in 2D case). Close to a gap node 
6*0 the order parameter is given by A{9) = A'{9 — 9q). 

We consider clean superconductors {A{9)t(9) ^ 1) so 
that the motion of quasiparticles inside and near the core 
is ballistic. In the quasiclassical approach the states are 



2 



described by the direction of the momentum 6 and the 
spatial position of the hne, on which the state is situ- 
ated. This spatial position can be described by the im- 
pact parameter b or, equivalently, by the angular mo- 
mentum I — kph. In d-wave voiiices bound states could 
extend far outside vortex cores .113 The characteristic ex- 
tension of these states in the direction of their momentum 
is of the order of the angle-dependent coherence length 
^{9) — VF/\^{d)\ where A(0) is the angle-dependent d- 
wave superconducting order parameter. 

Below the gap there is one chiral branch of bound 
states. Due to a special symmetry of (ID) BdG hamil- 
tonian in Andreev approximation for a state situated at 
a diameter {L^ 0) the eigenenergy of the bound state is 
exactly zero,E3 and the eigenfunction is 



'^{x) (X ^""^^ exp I — —J 



dyA{y,9) 



(4) 



with the extension of the order of (,{9). In Ref. the 
spectrum of states with impact parameters smaller than 
the vortex core size was calculated perturbatively, the 
correction to the energy being 



2/dx^exp(-^/dyA(y 



E = -b- 



/dxexp i-:^ JdyA{y 



(5) 



For small values of the impact parameter the logarithmi- 
cally divergent integral in the numerator of the previous 
expression should be cut off from below at the core size 
^0, that is, 



where the interlevel distance i£l 

uj„{e,o) = n„{e-9ofHi/\9-9o\) 



(6) 



(7) 

in a vicinity of a gap node ^^o — 2A''^ /vphp. Far 
from the nodes ujo{0) is of the order of JIq. However, 
for impact parameters larger than the core size the lower 
cutoff is given by the impact parameter itself. Therefore, 
to the logarithmic accuracy one finds that 



E{l,9)^-no{9 



9o)'ln 



(8) 



At large values of the impact parameter |&| 3> ^{9) the 
energy of the bound state saturates at the-gap value A{9), 
the asymptotic behavior being given by£3 



E{1,9) - -\A{9)\sign{l) + const^, 



(9) 



where const is of the order of unity. Saturation of the 
spectrum occurs at I ^ Ef/\A{9)\ ~ kp(^{9) where spa- 
tial variations of the order parameter along the quasiclas- 
sical trajectory are slow. 



In (charged) superconductors the effect of magnetic 
field on motion of quasiparticles leads to an additional 
term in the dispersion relation (||): SE — —luJ- This 
term is important for excitations with momenta close to 
gap nodes, where ^0(^*70) is small. 

Close to the nodes (k = ±a,±6) the extension £^(9) of 
the wave function of a localized state is very large and 
can exceed the intervortex distance Ry , which establishes 
the characteristic distance of the changes in the order pa- 
rameter and superfluid velocity in the vortex lattice. For 
angles 9 close to a gap node such that \9 — 9o\ < 9h, 
where 9h = A'/{vpRy) ^ £o/Rv ^ \/ H/ Hc2 , corre- 
sponding states are not localized any more and should 
be treated separately. For \9 — 9o\ ^ 9h the value ujo{9) 
becomes comparable to the cyclotron frequency Wc of the 
external field, which also contributes to uJo{9). We model 
this region by putting ujo{9) = luo{9h) ~ Wc for these an- 
gles close to gap nodes. The only characteristic of the 
spectrum in this region, which affects the value of the 
force, is the interlevel distance eo (see below). 

In addition to the just described branch of bound states 
there exist branches of the states with energies above the 
gap A(6'), that is, in the continuous spectrum. For a one- 
quantum vortex there is one chiral branch in the coittinu- 
ous spectrum, which crosses all energy levels once.E3 The 
interlevel spacing in the continuous spectrum is ujq = uJc- 

For later analysis we need to know the angle depen- 
dence of the interlevel spacing 



^,{9,E)^-[^ 



(10) 



for different energies E. For very small energies \E\ <ti 
A'9h this dependence is described by Eq.(^ above with 
the just discussed condition of saturation near nodes. For 
higher energies A'9h ^ \E\ < Aq the spectrum near the 
gap nodes is "continuous" : 



Ll!o{9,E) = LUc 



\9-9a\ < 9e 



(11) 



(where 9e is defined by |A(6l£;)| = \E\, 9e < 7r/4; 
9e = \E\/A' for \E\ < Aq). On approaching this re- 
gion from outside ujo{9,E) decreases. To investigate the 
behavior of luo{9,E) as we approach this vicinity from 
outside we use the expression (^ and find the relation- 
ship between I and 9 at given E as 9 ^ 9e- We find that 
I oc {9 — 9e)~^ 00 and luq ^ i}o{9 — 9e)^ ■ However, for 
I > kpRy the states are sensitive to the velocity fields of 
other vortices in the lattice, and therefore ujo{9, E) satu- 
rates at a finite value. This value is on the order of LOc^ 
which also contributes to dE/dl ([III). So, the behavior 
of ujo{9,E) as 9 approaches 9e from outside is quantita- 
tively similar to the behavior of ojo{9, 0) as 6* — > 0. 

As \E\ Aq the region of the continuous spectrum 
(O) grows and at \E\ > Aq the interlevel distance 
iLJo[9, E) ~ UJc for all 9. The difference in the behavior 
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of loq{9^ E) and, in particular, coexistence of the regions 
of continuous and discrete spectrum at the same energy 
i?, determines difference in the solutions of the KE at 
different energies (and therefore different temperatures). 

We can use the quasiclassical approximation as long 
as the temperature is larger than the interlevel spacing. 
The maximal value of L0a{9, E) is on the order of fio, i.e., 
we need T ^ f^o. 

Now, as we know the structure of the spectrum we can 
proceed to consideration of the dynamics of the distribu- 
tion function. 



III. KINETIC EQUATION FOR THE 
DISTRIBUTION FUNCTION 



We consider the problem in the frame of the vortex 
lattice, so that = 0. To calculate the effect of the 
SF (linear in the velocity of the normal component with 
respect to the vortex lattice v„ — v^) we put the exter- 
nal superflow = but keep the velocity of the normal 
component v„ finite. To the first order, the SF force is a 
linear function of v„ , posessing I?4-symmetry of the pair- 
ing state. Therefore, the force on a vortex is determined 
by two parameters: 



FsF = dgpCoK X V,; 



dgpCoKVn- 



(12) 



In particular, the longitudinal and transverse compo- 
nents of the force do not depend on the direction of v„. 
For simplicity we choose v„ = VnX. The dynamics of 
the distribution function (or the density matrix) can be 
described by the KE. Since 9 and I are canonically con- 
jugate, this equation reads (cf. Ref. 



dn{l,0) 
dt 



dE dn 
'de'dl 



dE dn 



(13) 



where n'^'' = np{E — v„k) = np{E{l,9) — Vnkp cos9) 
is the distribution function in equilibrium with the heat 
bath (the normal component), np is the Fermi distri- 
bution function. The treatment of the collisions in the 
relaxation time approximation is, of course, not accurate. 
Nevertheless, we believe that it qualitatively correctly de- 
scribes the physics of the system. The relaxation time r 
can be angle and energy-dependent, and we can easily 
incorporate this dependence into our results of this Sec- 
tion. However, for investigation of limiting cases we put 
t{9) — const. 

We study the stationary case and put dtn = in 
Eq. (|l3|) . If we change the variables 1,9 to E, 9, the equa- 
tion takes the form 



dn{E, ( 



n{E,9) - n'"'{E,t 



89 (^a{E,9)T 
The equation (|lj) can be easily solved: 



(14) 



n{E,9) = C{E)exp i J 



dip 



LUo{E, ip)T 



(15) 



[ — — nF(£'-w„fcFCOSx)exp I / 



dip 



ujq{E,lp)t 



The energy-dependent constant C{E) should be deter- 
mined from the condition of 27r-periodicity of n{E, 9) in 
9. 

Since we know the distribution function we can substi- 
tute it into the expression for the SF force (the force on 
vortex from the heat bath, which in the case of a crystal 
is the lattice with the impurities). Since the momentum 
of the excitations is given by:Q 

„ 1 [ d9dl , ^ . 

2 -^'^F[cos9;sin9)(n-no) 

[where nQ{l,9) = stands for normalization purposes 
and ensures that P = in equilibrium at T = 0; O(^) = 
for / < 0, 0(/) = 1 for / > 0; the prefactor 1/2 compen- 
sates for the double counting of particles |ajjid holes] , the 
contribution of the heat bath to the forceEj is fcf. ( 



' SF 



d9dl 



kp{cos9; sinl 



(16) 



To take into account all the branches we should sum up 
the expressions (^6|) for all of them. However, non-chiral 
branches never contribute to the force. 



IV. REACTIVE AND DISSIPATIVE FORCES IN 
VARIOUS LIMITS 

In principle one can substitute the solution (|l^) of the 
KE into the expression for the force (p^), and if the tem- 
perature dependence of the gap and the relaxation time 
is known, one obtains the temperature dependence of the 
SF force. However, it is interesting to investigate the be- 
havior of the force in different physical limits analytically. 

We are looking for an expression for the force to the 
first order in w„, so we expand the KE to this order. 
Namely, we characterize the deviation of the distribution 
function from its equilibrium value at w„ = by a func- 
tion I/: 



i{E,9) ^ npiE) 



dnp 
~dE 



Vnkpi^iE, 9). 



The stationary KE takes the form 



dt 



89 LUoiE,9)T 



(17) 



(18) 



where v'^'^ — cos 9. This equation has a unique solution, 
which is an odd function of 9: v{9 -f- tt) = —v{9). The 
quasiclassical equations of motion 
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dE 



I 



dE 
~d9 



(19) 



describe the motion along a curve of fixed energy in (6, l)- 
plane, and the KE ( [T8| ) keeps track of variations of the 
distribution function v along this curve in the stationary 
case^ The time of motion along this trajectory is given 



dip 



ujo{E,ip) 



(20) 



and the KE describes the relaxation of v towards its equi- 
librium (^-dependent) value: 



'dt 



eg 



(21) 



The solution of this equation depends on the function 
i/^'^{t) and the relaxation time t. 

The total time of motion along the closed trajectory 



to{E) = / 
Jo 



dip 



uJoiE.if) 

determines the structure of the ener|g¥|ievels. The quan- 
tization of the quasiclassical motiorJlj'O shows that in a 
vicinity of energy E the spectrum is equidistant with the 
interlevel distance 



eoiE) 



27r 



(22) 



For small \E\ ^ ^'Oh the interlevel spacing equals 
eo = eo{E = 0) and is given by Eq.(||). 

The solution ( p^ can be simplified in various limits. 
The behavior of the force is different in three regimes: 
moderately clean with w™°^t ^ 1, "very clean" with 
eoT ^ 1 <C and "extremely clean" when cqt ^ 1. 
Usually the term "superclean" is used for the regime with 
3> 1. However in the mixed state of d-wave ma- 
terials vortex dynamics is different in two subregimes. 
Therefore we use the terms "very clean" and "extremely 
clean" for these two subregimes. Note that egr depends 
on external field H and temperature, while Qqt is de- 
termined by temperature only. In subsections A, B, 
C below we consider these three regimes at low tem- 
peratures T ^ TcOh [this condition is consistent with 
T » fio if H/Hc2 > (Ao/E'f)^]. At corresponding en- 
ergies <C A'S/f the quasiclassical interlevel spacing is 
given by ujq{6,0). The case of higher T is discussed in 
subsection D. 



A. Moderately clean limit 



In the moderately clean case the parameter ujqt is 
small for all values of 9: w^^^t ^ 1, and v^'^{t) is a slow 
function of t (on the scale of t). Therefore, the solution 



V of the KE everywhere is very close to v'^'^ . Simple cal- 
culation shows that the value of the SF force in this limit 
coincides with that in the s-wave case. The contribution 
of the bound states is given by: 



P^bound _ 1 [ dOdl . . 



2kCqv„ 



dEde 
2n 



dnp \ 
"dEJ 

/dO 
— (008 6*; sin 6*) sin^ 

Anp{-\m\)-npme)\)]. 



cost/; smW smi 



(23) 



In the continuous spectrum only the chiral branch con- 
tributes: 

/do 
— (cos 6*; sin 6*) sin 61 
27r 

^[l-np{-\m\)]- (24) 
So, the total force is 



' SF 



Cqk X v„ 



4f = 1- 



(25) 



The first correction (in small t) gives the dissipative 
force: 



dlp^T{u;oie,0)), 



(26) 



where the angle brackets denote averaging over angles 9. 

In this limit the SF force almost exactly cancels the 
Magnus force, therefore the total non-dissipative force 
(and the Hall conductivity) is given by the force, propor- 
tional to the difference p ~ Cq (see Introduction) and a 
correction to fE5h: 



5d^ 



SF 



(27) 



So, the total dissipative and non-disspative forces are of 
the same order as in an s-wave superconductor with the 
same value of the gap Aq. 

The result (^5|) is independent of the structure of the 
vortex core. However, the corrections (p^), ( p7| ) depend 
on this structure. 



B. Very clean limit 

In this limit eqt <C 1 ^ t^p'^^r the value of ujqt is 
small near the nodes but it is large far from the nodes. 
t{0) is almost constant far from the nodes, but it varies 
rapidly (on t time scale) in their 0^-vicinities where 6^ = 
{^Iot)~^. In an internode region 9o + a < 9 < 9q + tt /2 — a 
(where 9r <^ a <^ 1 and 6*0 = 0, 7r/2, tt, 37r/2 is a gap 
node) the point in (0,l)-plane, describing the quasiclas- 
sical motion (pj|), spends an amount of time small com- 
pared to the relaxation time r. It follows from the KE 
that v is almost constant in such a region (variations of 
V are much less than unity; we assume, and this assump- 
tion is justified by the form of the solution below, that 
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everywhere on on 0-circle v is of the order of 1 or less). 
On the contrary, as the point in (0; /)-plane approaches a 
node its motion gets slower and v adjusts to a local value 
of v^'' (the positive i-direction is the negative 6'-direction, 
therefore the jump in the value of v occurs on the right- 
hand side of a node, see below). In a small 0- vicinity of 
a gap node v'^'' is almost constant and the solution of the 
KE (21) can be representectil as y = A + Be~*/'^ with 
constant A and B. In the very clean limit, the change 
in the exponent across a gap node to/4T — 7r/2eoT is 
very large, and matching the solutions in two different 
regions leads to the conclusion that in the internode re- 
gion [^o; 6'o+7r/2] one has v = z/'^'?(6'o +7r/2), and v jumps 
between two constant values near the nodes. In terms of 
V the force is given by 



¥sF = 2kCow„ 



dEdO 
27r 



dnp \ 



X (cos^?; sin t 



UJQT 



(28) 



Using the KE we can substitute the last fraction in the 
rhs by dv/dO. This angle derivative of v has sharp peaks 
near the nodes, and simple integration shows that in the 
very clean regime 



d^ -d'l 



(29) 



C. Extremely clean limit 

In the extremely clean limit r is very large, namely, 
tQT ^1. It means that the period of a quasiclassical 
trajectory Iq is much less than r. Therefore, v does not 
have enough time to change its value during the motion, 
i.e., V = const. Since the solution of the KE ( ^l|) should 
be an odd function of 0, this constant is in fact 0. So, to 
the first order the total force is given by 



FsF = 2kCow„ 



dEde 
2tt 



dn 



—— (cos 6*; sin 6*). 

oE J uqt 



At low T ^ Tc9h the leading contribution to the force 
in this limit is dissipative: 



(30) 



eo is the only characteristic of the spectrum, which enters 
the expression for the force in this limit. 



D. Higher temperatures 

Moderately clean limit. At higher temperatures 
T > TcOh one should account for energy dependence of 
ujo{9,E) at characteristic \E\ ^ T. In the moderately 



clean regime a;™°^r is still very small (this parameter 
depends on T only through Ag and r) and the force is 
given by Eg. (p5|) . To obtain corrections one should sub- 
stitute ujo{0,(^m Eqs.(|2|),(||) by lJ(){0,E) and average 
over energies.E2l At low T <$^Tc these corrections remain 
of the same order of magnitude as at T ^ TcOh- As T 
approaches Tc they decrease down to lOcT and —{uJc't)'^, 
respectively. 

Super clean limit. In the superclean regime 
^max^ ^ ^ ^j^g situation is more complicated. The time 
of quasiclassical motion along the closed trajectory at a 
given energy E ^ A.'9h has two contributions. The time 
of motion outside ^s-vicinities of gap nodes is on the or- 
der of 27r/eo (cf. Section 2), while the motion inside these 
vicinities takes time on order of We/^c 3> 2-K/tQ. 

Low fields. If e^T <^ 1 then the sample is in the very 
clean limit (subsection B) at all temperatures, i.e., the 
system ( [ig| ) spends a large amount of time (in units of r) 
near the nodes and a small amount of time between the 
nodes. Therefore, the solution of the KE is of the same 
form as in the subsection B. In 9e vicinities of the nodes 
ly = v^'i. Using the KE ( pi] ) it is easy to show that on the 
interval [6q + Oe'tOq + 7r/2 — 9e\ the solution is given by 
V = v'^''{9q -|-7r/2 — Oe) to the accuracy o(l). The transi- 
tion between these two types of behavior of the solution 
occurs in an intermediate region of the width 6^ 

There are two contributions to the integral over 9 in 
(^): one from ^^-vicinities of the nodes, and the other 
from the boundaries of these regions where u jumps to 
its equilibrium value in the internode region. The inte- 
gration shows that (to the accuracy o(l)) 
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4-tanh^. 

2T 




(32) 



The last line of Eq. ( p2[ ) represents the contribution of en- 
ergies > Aq and is exactly the same as in the s-wave 
case (cf. Ref. |). 

At very low temperatures T <C Aq these expres- 
sions give ([29|). We see that the SF force is almost T- 
independent in this region of temperatures. In the oppo- 
site limit T > Ao(T) (which corresponds to \Tc - T] < 

Tc) we have = and dgp = 1 regardless of the value 
of uj]J^"'^t: in this region the SF force exactly compen- 
sates for the Magnus force, and the vortex does not ex- 
perience any force in the limit of the normal state (note, 
however, that our approach is not applicable to a very 
narrow vicinity of Tc where Aqt ^1). 

Higher fields. In the other case cqt 3> 1 the time 
of periodic motion is small at sufficiently low temper- 
atures T <^ TcLOcT. In this limit the situation is sim- 
ilar to that of subsection C, and spectral flow is sup- 
pressed. To obtain an expression for the force one should 
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average over energieaiS Eg. (pOD with ep substituted by 
eo{E) = ttlOc/^Oe- Evaluation shows that: 



(38) 



''SF 



T 1 

Tr LOrT ' 



T 



TrUJrT 



(33) 



At low T <^Tc one can calculate prefactors in these ex- 
pressions: 



where Pmn are matrix elements of the operator of linear 
momentum. 

At low energies E = —luqIO)! and quantization leads 
to a hermitian operator H = iujQ{9)de + ideojo{0) /2 as a 
hamiltonian. The eigenfunctions of Ti. are 



„ _81n2T 1 



If LUcT 3> 1 the temperature in this extremely clean 
regime can reach the values T ^ Tc- At these tempera- 
tures the contribution of the states above the gap Aq to 
Sgp,dgp is of the same order as the contribution of the 



states below the gap. At T ~ Tc expressions (^3|) corre- 
spond to conductivities of an s-wave superconductor at 
T ~ Tc or a normal metal in the limit uJcT ^ 1. 

If ujcT ^ 1 ^ eo''' another regime is possible: at higher 
temperatures TcOJcT <^ T < Tc the quasiclassical motion 
near the nodes becomes slow (on t time scale), the sys- 
tem is in the very clean regime, and the force is described 
by Eqs.dsTl), (|32 



V. ON THE DENSITY MATRIX 

The parameter cqt is determined by the interlevel dis- 
tance, while the relation of Qqt to the structure of the 
spectrum is unclear. In this Section we consider the gen- 
eral KE for the density matrix in terms of exact (not qua- 
siclassical) eigenstates of the hamiltonian. We demon- 
strate a solution of this equation and evaluate character- 
istic parameters determining the strength of the SF. We 
discuss the origin of the parameters floT and e^r. 

The KE for the density matrix p is 



dp 
dt 



i[p,n] 



fig 



(34) 



In the stationary case 



rrnn 



and therefore 
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Er, 



(35) 



Here n, m enumerate the eigenstates of Ti.. In two ex- 
treme cases we have: 



-p^'Jr, if \E, 
The SF force on the vortex is 



if |^™-^„|t«1, (36) 
£;„|r>l. (37) 



fNiO) = 



eo 



M0) 



exp 



^iE 



N 



dip 
uj^iip) 



(39) 



(cf. Refs. where a non-hermitian operator iuji^de 

was used), Eff = NtQ are determined by periodicity of 
the wave functions. Matrix elements of the momentum 
operator P„m [which also determines the density ma- 
trix p'^'' = np{'H — v„P)] depend only on the difference 
N ^ m — n. In terms of P+ — Px + iPy the force is given 
by 



SF 



N 



iN 



1 + iNe^T 



(40) 



at low T ^ Tc\/ H / Hc2 (one can split summations over 
n and -^V— j= m — n in ( |3q ) provided that T ^ fio)- To 
evaluateE^ the force at higher temperatures one should 
substitute eo by eo(i? ^ T) ^ lOcTc/T. Matrix elements 
of P+ are 



P, 



N 



dx 



where 6{x — —eat) is defined by Eq.(pO|), and are non- 
zero only for iV = 4fc + 1. Analysis shows that for small 
TV (such that - 1) one has P+ = 2V2{-l)''kp /ttN . 
P^ deviates from this expression at N ^ f^o/eo ~ 
(these N are on the order of the maximal value of dO/dx) 
where it slightly increases, and P^ rapidly decreases for 
\N\ ^ 9^. So, this dependence determines two energy 
scales: the interlevel distance eo and eo(Oo/eo) = f^o- 
Analyzing Eq.(pO[) one can single out the same three 
regimes as in the previous Section. 

In the moderately clean regime one can neglect the sec- 
ond term in the denominator in (40). Two regions of N 
contribute to (only reactive) force: the region of small 
\N\, which gives a contribution 2/7r to dgp, and the re- 
gion A'' ~ /eo contributing 1 — 2/7r to dgp . 

In the very clean regime eoT <C 1 <C ^qt the second 
term in the denominator of Eq.(^) suppresses the con- 
tribution of the region N ^ ^l^ I as this term is on the 
order of iloT S> 1 in this region. On the other hand, 
N ~ (eor)^^ contribute to the dissipative force in this 
regime, the contribution to Sgp being 2/7r. Finally, in 
the extremely clean regime egr ^ 1 the denominator in 
( ^ ) suppresses all the terms, and to the first order the 
force is Fsi? — P'^'^/r. Calculation of the equihbrium 
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value of the linear momentum (at given v„) leads to the 
result @). 

In the s-wave case 7^ only for = 1, and Eq.pO|) 
gives the same result as in Ref. ^. 

Thus, we found that the reactive force is due to transi- 
tions between the levels separated by distances a) on the 
order of interlevel distance, and b) on the order of max- 
imal value of ujq{0). However at finite relaxation rate 
1/t only the terms with energy difference Neo < l/r in 
( [40| ) are effective. The dissipative force, naturally, comes 
from transitions over c) the distance l/r. Contributions 
b) and c) compete with each other: the former is sup- 
pressed if fipT 1, and the latter in the opposite hmit. 
In the extremely clean limit all the transitions are sup- 
pressed. 



VI. SUMMARY AND DISCUSSION 

Our analysis shows that like in s-superconductors, in 
d-materials the SF along anomalous branches of the spec- 
trum of fermions bound to vortex cores produces a contri- 
bution to the force on a moving vortex. In the regime of 
extreme SF (fast relaxation in the core) the value of the 
non-dissipative force is the same as in the s-wave case. 
This result could be anticipated: ki fact, consideration 
of the action for a superfluid shows£l that the SF force is 
proportional to the number of the states below the Fermi 
level (Co / m) regardless of the pairing symmetry. 

In the superclean regime coexistence of bound states 
with states in the continuous spectrum leads to new fea- 
tures in vortex dynamics. As the relaxation time grows at 
constant T and H (the system gets cleaner) , the spectral 
flow contribution to the reactive force is gradually sup- 
pressed. The motion of a vortex with respect to the nor- 
mal component (crystall lattice) produces a perturbation 
of the Bogolyubov-de Gennes hamiltonian vP leading to 
transitions between energy levels. Transitions between 
the levels n and m do not lead to a net loss of momen- 
tum of quasiparticles if the time derivative En — E„i of 
the phase of an off-diagonal element of the density ma- 
trix pmn is much larger than the relaxation rate l/r. 
Such transitions effectively contribute to the force only 
for \En — Era\T < 1. In s-wave vortices the operator of 
linear momentum leads to transitions only between the 
neighboring bound levels with energy spacing wq, and 
the parameter ujqt determines the total force on a vortex 
(cf. Ref. ^ . In the d-wave case the perturbation of the 
hamiltonian effectively leads to transitions between levels 
with energy differences in the range [eq; ^o\- So, features 
of vortex dynamics depend on two parameters, Qqt and 
eoT. The reactive force is completely suppressed only in 
the extremely clean limit. 

The dissipative contribution is small in the limits of 
very large (extremely clean limit) and very small (mod- 
erately clean limit) r and it acquires its maximal value in 
the intermediate (very clean) region. In this limit dissipa- 



tive and reactive contributions are equal and independent 
on r and the details of the spectrum. 

In connection with this discussion of different symme- 
tries of the pairing state let us mention that our calcu- 
lations are not sensitive to the sign of the order parame- 
ter. Only the anisotropy of the absolute value of the gap 
counts. So, the results are equally applicable to d-wave 
and extremely anisotropic s-wave superconductors (with 
gap nodes). 

Let us discuss some conditions and restrictions we used 
in our calculations. We have mentioned above that sim- 
ple quasiclassical description of Andreev bound states is 
valid in the clean limit A{6)t ^ 1. In the moderately 
clean limit the main contribution to the force is due to the 
states away from the nodes. Therefore, the contribution 
of vicinities of the nodes is inessential, and the condition 
Aor 3> 1 is sufficient. In the superclean limit the main 
contribution comes from vicinities of the gap nodes of 
the width max{Or,dE)- The condition A{9t)t ^ 1 can 
be translated to EpT 3> 1, which is even weaker than 
Aor > 1. 

In our calculations we omitted the logarithmic factors 
[cf. Eqs.(^,(^)]. In fact these factors are unimportant 
in the moderately clean limit, where mostly the regions 
away from the nodes contribute, and in these regions the 
logarithms are absent (see Section 2). In the superclean 
limit, where the contribution of the states close to the 
nodes is essential, the logarithmic factors could only af- 
fect the boundaries between the different regimes. It fol- 
lows from the derivation in Sections 4,5 that they could 
not change the results in the very clean regime. As for 
the extremely clean limit, the logarithms enter the result 
only through eo- 

We studied 2D case so far. To account for the third 
dimension (with H |j z) we can integrate over all modes 
of motion along the third axis z. z-dependence of the 
spectrum of the bound states is slow (within an order 
of magnitude), therefore our results would not change 
qualitatively in 3D. 

Using our results one can investigate temperature, 
magnetic field and purity dependences of transverse an 
and longitudinal cfq conductivities, which are given by 
(cf. (0) and Ref. ||): 



,|| 

— rP" 



SF 



mH 



(41) 
(42) 



In the moderately clean limit the behavior of dissipa- 
tive and non-dissipative forces is analogous to that in 
s-superconductors. In the superclean limit the values of 
the Ohmic and Hall conductivities and their relaxation 
time, magnetic field and temperature dependences differ 
from those in the s-wave case (at the same Ao). 

In the moderately clean regime d}gp and 1 — dgp are 
small and almost field independent, aH,ao oc H~^, 
an *C (tq and the Hall angle is very small. The su- 
perconductor is always in this regime {A^t/Ef ^ 1 or 



8 



Tc — T<^ Tc) just below Tc, and this region could extend 
down to some finite or zero [if A§(T = 0)t{T = 0)/Ef < 
1] temperature. In the latter case the superclean limit is 
unreachable. 

In the very clean limit the conductivities are given by 
Eqs.{|l|),(i|) and Eqs.(|ll),(||). This regime is a sub- 
regime of the superclean one with an additional condition 
that eoT < max(l; y/Hc2/HT/Tc). If in addition T < 
then ao and cth are T and r independent and cx H^^ 
[cf. Eq.(p9|)]. Experimentally the superclean regime in 



high-Tc compounds was reached, e.g., in Ref. 25 



In higher fields [eoT > max{l: Hc2/ HT/Tc); ex- 
tremely clean limit] ctq is small and the Hall angle is close 
to it/2. In this regime ao oc 77~^max(T/Tc; ^/WJlI^). 
In all the three regimes the conductivities depend on T 
directly and through Aq and r (see Section 4) . The exact 
view of the phase diagram in T-H plane is determined 
by parameters of a sample, in particular, by temperature 
dependence of the relaxation time r. 

It would be interesting to investigate the influence of 
anisotropy of the Fermi surface on vortex dynamics. In 
the moderately clean limit ajiiaotropy should not affect 
mutual friction parameters. Q'E3 However in the super- 
clean regime it may be important. 

In conclusion, we presented an analysis of the contribu- 
tion of the quasiparticles bound to vortex cores to vortex 
dynamics in the mixed state of d-superconductors with 
lines of gap nodes. We found three regimes of vortex 
dynamics. The behavior of Hall and dissipative (Ohmic) 
conductivities in the moderately clean regime is similar to 
that in s-wave superconductors, while in the superclean 
regime new features were found. 

As this article was in preparation, I learned about a re- 
cent work by N.Kopnin and G.Volovikcj devoted to the 
same subject. We note that our results for finite temper- 
atures T > Tc^H/Hc2 differ from those of Ref. H. 
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